Measuring equality horizon with the zero-crossing of the galaxy correlation function 
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The size of the horizon at the matter-radiation equality is a key scale of the Big Bang cosmology 
that is directly related to the energy-matter content of the Universe. In this letter, we argue that 
this scale can be accurately measured from the observed clustering of galaxies in new large scale 
surveys. We demonstrate that the zero-crossing, r^, of the 2-point galaxy correlation function is 
closely related to the horizon size at matter-radiation equality for a large variety of flat ACDM 
models. Using large-volume cosmological simulations, we also show that the pristine zero-crossing 
is unaltered by non-linear evolution of density fluctuations, redshift distortions and galaxy biases. 
This makes rc a very powerful standard ruler that can be accurately measured, at a percent level, 
in upcoming experiments that will collect redshifts of millions of galaxies and quasars. 
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PACS numbers: 98.80.-k, 98.80.Es,98.65.Dx 

The science exploitation of upcoming experiments such 
as BOSS, DES, LSST, BigBOSS and Euclid that will sur- 
vey the uncharted universe taking spectra and images of 
millions of galaxies and quasars, combined with the anal- 
ysis of Grand Challenge cosmological simulations and de- 
velopment of theoretical models, will be critical to accu- 
rately measure the properties of our Universe (e.g. [l| for 
a review). 

The physics governing the evolution of perturbations in 
CDM universes imprints two distinct length scales in the 
galaxy distribution which are useful as "standard rulers" 
for distance estimates: (1) the scale of the particle hori- 
zon at matter-radiation equality and (2) the sound hori- 
zon scale at the baryon-drag epoch before recombination 
[3, [^. The Baryonic Acoustic Oscillations (BAO), re- 
lated with the later, have been detected in the SDSS, 
2dFGRS and WiggleZ galaxy power spectra at a char- 
acteristic scale khao > 0.05/iMpc~^, and as a single 
peak in the galaxy two-point correlation function (CF) 
at rbao « 105/i-iMpc (e.g. 

The total physical matter density (wm = ^mh^ — 
flcdmh^ + ^hh^) and, thus, the moment of equality can 
be measured from CMB fluctuations alone, but the mea- 
surement errors decline almost by a factor of two if BAO 
clustering estimates are included @. It would be helpful 
to have another method for testing the equality epoch, 
which is more sensitive to than BAO. 

The particle horizon at matter-radiation equality is 
associated with the most prominent feature of the lin- 
ear power spectrum of density fluctuations P{k) - the 
turnover or maximum - at the characteristic scale fcmax 
~ 0.015 ft.Mpc""'^. This corresponds to the transition from 
P{k) ~ k for a Harrison-Zeldovich scale-invariant spec- 
trum to a P{k) ~ k~^ spectrum due to modes that were 
not growing after they entered the horizon during the 
radiation dominated era [l^. Thus, we expect that the 
location of this turnover will be related to the scale of 
the horizon at matter-radiation equality k^q Q. 



The equality wavenumber k^q is given by k^q = 
) with Ooq = {ujj/ujn^){l + 0.2271 A'"cff) and 
being the scale factor 
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H, 

and expansion rate at matter-radiation equality 2]. We 
adopt a radiation density = Vl^h"^ = 2.469 x 10~^ and 
iVoff = 3.04 for standard neutrino species jSl]. Hence, it 
can be seen that the equality horizon scale depends solely 
on the matter density The size of the hori- 

zon at the matter-radiation equality epoch is defined by 
the comoving distance roq = I?h ~ jl^°'*[a^'^ / H{a)\da — 
(4 — 2\/2)k~^. Thus, roq scales with matter density 
as Toq cx . We obtain fccq ~ 0.014 /iMpc~^ and 
Tcq ~ 85ft.~^Mpc adopting the latest matter density ob- 
servational constraints 0. 

The turnover in P{k) shifts to high wavenumbers for 
models with larger matter density maiidy as a result of 
the scaling relation between the equality wavenumber k^q 
and Win- The matter perturbations remain frozen after 
entering the horizon, but since this process is not instan- 
taneous the transition from the P{k) (x k to P{k) oc k~^ 
is broad over about an order of magnitude in k. For dark 
matter-baryonic models the expected scaling with matter 
density might involve an additional dependence on other 
cosmological parameters. In order to derive a relation 
between fcmax at the maximum of P{k) and the equality 
scale fccq, we computed the linear matter power spec- 
trum of density fluctuations of a flat ACDM universe, 
using CAMB 11|, for a number of cosmological mod- 
els with varying and around our flducial model. 
For this fiducial model, which is consistent with WMAP7 
cosmological parameters and other recent cosmological 
constraints (i93 and references therein), we assume the 
matter density = 0.132, the normalization of fluctu- 
ations (78 = 0.82, the primordial index of scalar pertur- 
bations rig = 0.95, the baryon density Wb — 0.023, the 
Hubble constant Hq = lOO/ikms'^Mpc-i with h = 0.7, 
and iVcff = 3.04. In Fig. [1] (top panel), we show the lin- 
ear matter power spectrum of our ACDM models, com- 
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FIG. 1. Linear power spectra P{k) (top panel), and their 
corresponding correlation functions f(r) (bottom panel) for 
ACDM models with a fixed baryon density Uh = 0.023. Thick 
lines corresponds to our fiducial model. When ujm increases, 
the turnover in P{k) and its characteristic scale fcmax moves to 
larger k. Thus, the zero-crossing in the correlation function 
^(rc) = 0, located beyond the BAO peak, shifts to smaller 
radii. For each model, both P{k) and ^(r) have been normal- 
ized to its maximum and BAO peak respectively. 

puted with CAMB, for different uj^ values and kept fixed 
the other cosmological parameters that determine P{k), 
to that given by our fiducial model. As expected, the 
turnover location fc^ax is close to fceq and scales with Wni- 
Yet, the precise turnover scale fcmax depends slightly on 
other cosmological parameters such as the primordial tilt 

Us. 

We show in Fig. [2] (left panel) the dependence of fcmax 
as a function of the expected equality wavenumber fcoq for 
different models with fixed baryon density 25% below and 
above that of our fiducial model with Wb = 0.023 (solid 
line), i.e. ^ 6a around best observational constraints. 
Matter density a;,„ varies over a large range from 0.05 to 
0.25 around our fiducial value Wm = 0.132 (solid symbol). 
We kept fixed as and Us- The dotted line represents the 
fcmax = fcoq relation which scales only with Wm as fcoq = 
0.104ajin /iMpc"^. The difference between fcmax and fcoq 
is less than 20% for the range of parameters considered 
here. Models with smaller baryon content show a better 
agreement between fcmax and fceq. The weak dependence 
of fcmax (fcoq) On baryon density lu^ can be approximated 
by 

_ / 0.194 \ , 0.685-0. 121 1ogio(<^b) (i\ 
/i-max — ^0.321 J '^eq ■ \'^) 



We also analyzed models with different primordial tilt rig 
^ 2.5ct around the latest constraints 9], and found only 
a small dependence of fcmax on Ug. 
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FIG. 2. P{k) turnover location fcmax as a function of the 
horizon scale at matter-radiation equality fceq (left panel), and 
its relation with the corresponding zero-crossing position Vc 
in the linear ^(r) (right panel) for ACDM models with three 
different baryon densities. Solid symbols are for our fiducial 
model. 

The existence of the turnover in the galaxy power spec- 
trum has been discussed extensively in the literature, al- 
though its detection has not been reported up to now 
in any completed galaxy redshift survey such as the 2dF- 
GRS, SDSS or WiggleZ Yet, little attention has 

been paid to the imprint of the horizon scale at matter- 
radiation equality on the zero-crossing in the galaxy two- 
point correlation function ^(r). Although the very small 
clustering signal at those large scales is severely affected 
by potential sources of systematic errors [ll,[l|, here, m 
this letter, we highlight that the zero-crossing deserves 
momentous attention both from theory and observations. 
In the early 90s, Klypin & Rhee [l^ proposed to use the 
zero-crossing measurements of the galaxy cluster CP as a 
sensitive test for the shape of the power spectrum of ini- 
tial fiuctuations. See also [l^ for a detailed discussion on 
how, in their case, a sharp maximum in the power spec- 
trum would be mapped into the correlation function. 

The detection of the zero-crossing in ^(r) is a funda- 
mental prediction of ACDM or any other dark matter- 
baryon cosmological model [13, [U [3 ■ The correlation 
function is related to P{k) by the Fourier transform p^ : 

The correlation function must have a zero-crossing 
iifc) = at some radius because P{k) — )• when fc — )■ 0, 
implying that £,{r)r^dr — 0. Note that for our fidu- 
cial model, the zero-crossing Vc happens only at a scale 
^ 20% larger than the position of the BAO peak, i.e. at 
129.3 h-^Mpc (see Fig.[Tl bottom panel). 

However, the linear ^(r) for models with Wm ^ 
0.2 shows a first zero-crossing located at scales Tc ^ 
SO/i-^Mpc, i.e. below the BAO feature (see Fig. [J). 
We recaU that the SDSS and WiggleZ galaxy CFs dis- 
play a positive clustering signal well above ~ 80 /i^^Mpc, 
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FIG. 3. Top curves show the zero-crossing scale rc as a func- 
tion of matter density cj~^ and equaUty horizon size r^q (top 
axis) for our ACDM CAMB models with fixed values of cjb- 
This relation can be well represented by a broken power-law 
as given in Eq. ((3|. As a reference, we also show the acoustic 
peak position rbao (bottom curves). Solid symbol corresponds 
to our fiducial model. Zero-crossing estimates and 1-cr errors 
obtained from best-fit and a stretch parameters from dif- 
ferent survey CF data are shown (open symbols). The inset 
panel shows the corresponding observational rcq/Dviz) esti- 
mates for different redshifts. Also displayed are predictions 
for ACDM models with ojm =0.11 (dashed line), 0.132 (fidu- 
cial, thick solid line) and 0.17 (dotted line). 



at least up to the BAO peak located at ~ 105ft.~^Mpc, 
which aheady imposes a matter density upper limit of 
Wm ~ 0.2 (the stacked ga laxy CF of these surveys gives 
a BAO detection at 4.9(7 [iMl)- Here, we study the zero- 
crossing in ^ (r) located beyond the BAO feature for mod- 
els with 50% matter density variation around our fiducial 
value, i.e. 0.066 < < 0.198. 

Fig. [5] (right panel) shows the relation between the 
zero-crossing position Tc in the linear regime and its cor- 
responding P{k) turnover location at fcmax for models 
with 0.066 < Wm < 0.198. It is clearly seen in Fig. [3]that 
Tc is determined by the scale of the horizon at matter- 
radiation equality. Note, that Tc shows a much steeper 
increase with oj^^ as compared with the position of the 
BAO peak rbao- Moreover, the BAO peak position shows 
a larger dependence on baryon density Wb, e.g. about 7% 
and 3.5% variation in rbao and Tc respectively at fiducial 
Wm = 0.132 (solid point). The zero-crossing displays a 
small dependence on Wb as can be seen in Eq. [TJ For 
our grid of CAMB models with 0.066 < w,n < 0.198 and 
0.0172 < Wb < 0.0284 we obtain the following relation 



between the zero-crossing and the comoving size of the 
matter-radiation equality horizon, which is well described 
by a broken power-law, i.e. 



''oq.O 



(3) 



where r, 



eq,0 



129.5 /i-^Mpc, bi = 0.425, d = 6.5. Pa- 



rameters A = 14.4780;^"" '''"'' , and 62 = 34.456 Wb+0. 171 
encode a slight dependence on baryon density. The size 
of the equality horizon r^q depends only on the equality 
characteristic scale fcoq through the matter density, i.e. 
Teq = (4 - 2%/2) k-^ = 11.231 cj^i h-^Mpc. 

To estimate the value of Tc with current galaxy surveys 
we use the constraints on the stretch parameter a(zoff) = 
Dvizcs)/Dv^f,dizeB) (Dviz) = [(1 + zfD\cz/Hiz)]^/^ 
is the dilation distance [4] and Da is the angular diam- 
eter distance), obtained by different groups from their 
best-fit cosmology modeling to the observed SDSS and 
WiggleZ redshift-space ^(s). Typically, this model fit 
is done on the redshift-scale range 20/i^^Mpc ^ s ^ 
180/i~^Mpc. This allows us to estimate the expected 
position of the zero-crossing using the relation rc = 
ct{zcs) fcfid- -Dv,fid(z) and r^^d are the dilation distance 
at Zeff (survey mean redshift distribution) and the zero- 
crossing scale for the fiducial cosmology. Fig. [3] shows 
our estimates taking constraints on a( zp,ff) and uj^ from 
SDSS-LRG E05 3: SDSS-LRG DR7 Ullll, WiggleZ 
& Stacked WiggleZ+SDSS+6dFGS [l|; and SDSS-LRG 
DR6 combined with CMB 0. 

The accurate modeling and interpretation of both P{k) 
and ^(r) galaxy clustering statistics present some diffi- 
culties. The shape of the linear matter power spectrum 
P{k) is distorted by the nonlinear evolution of density 
fluctuations, redshift distortions and galaxy bias even at 
large-scales k < 0.2/iMpc~^ [2l| - [2i |. In this regard we 
should also draw attention to the impact of those in- 
evitable effects on the expected zero-crossing scale in the 
redshift-space galaxy correlation function, ^(s), as im- 
printed by the horizon scale at matter-radiation equality. 

We show in Fig. 2] redshift-space ^(s) at z = 0.5 (close 
to the mean redshift of the BOSS and WiggleZ galaxy 
distributions [8, 25]), drawn from our MultiDark and 
LittleMuk suite of A^-body simulations of a flat ACDM 
cosmology with the same parameters as for our fiducial 
cosmological model. The MultiDark simulation (MDRl) 
was done using the ART code. It has 2048^ particles 
in a l/i~^Gpc box (see [2^ and www.multidark.org for 
details). The mass and force resolutions are 8.72 x 
IO^/i-^Mq and 7ft."^kpc. LittleMuk (LMuk) consists 
of three GADGET [H realizations of 1280^ particles of 
a larger box with 2.5 h~^Gpc on a side. Force resolution 
and particle mass are 28 h~^kpc and 5.58 x 10^^ /i~^Mq. 
The LMuk total volume 46.875 {h^^Gpcf allows us to 
compute the mean and the variance for the estimates of 
£_{r) for dark matter halos. Dark matter halos (and sub- 
halos) were identified with the Bound-Density-Maxima 
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FIG. 4. Redshift-space correlation functions at z = 0.5 for 
three halo ("galaxy") samples with different bias in the Mul- 
tiDark and LittleMuk ACDM simulations. Results are com- 
pared with the linear ^(r) (thick solid line) that corresponds 
to our fiducial cosmological parameters. The inset panel 
shows the real-space correlation function. The zero-crossing 
in real space is the same regardless of galaxy bias. There 
is a small ~ 0.4% shift in the redshift-space. Dotted lines 
are statistical errors estimated for "galaxies" in the volume of 
46.875 (/i-^Gpc)^ 



(BDM) algorithm |26|]. We then use a simple, non- 
parametric abundance-matching prescription, to connect 
dark matter halos (and subhalos) in our simulations with 
galaxies by selecting them above a given maximum cir- 
cular velocity Knax- This procedure is able to predict the 
clustering properties, and therefore the two-point CF and 
halo occupation distribution of observed galaxies for dif- 
ferent number densities (e.g. [lEIlll)- Knax thresholds 
for the three "galaxy" samples with increasing biases in- 
dicated in the plot are 180, 350 and 600 kms~^ roughly 
corresponding to the Emission Line Galaxies (ELGs, 
b = 1.4), Luminous Red Galaxies (LRGs, b = 2.0) and 
QSOs (6 = 3.0) in the major experiments discussed here. 

Results presented in Fig. |4] evidence that the zero- 
crossing position Tc in the dark matter halo CF appears 
at the location predicted by linear theory. Thus, Tc is 
not distorted by the non-linear evolution of density fluc- 
tuations, redshift distortions and bias. This makes the 
zero-crossing an invaluable tool to constrain the size of 
the horizon at matter-radiation equality and the shape of 
the primordial power spectrum of initial fluctuations. We 



measure Tc = 129.8l5'g /i~^Mpc from the mean LMuk 
"galaxy" sample (thin solid line in Fig. |4]) with errors 
scaled according to statistics similar to that of the SDSS 
LRGs (^ 10^ galaxies). Zero-crossing uncertainty goes 
down to a small value of ~ 1.5% if we consider a sur- 
vey with 10 times more galaxy statistics such as BOSS 
(dotted lines) or ~ 0.3% in the case of 20 million galax- 
ies, as it will be targeted by the upcoming new surveys 
like BigBOSS or Euclid. Our estimates are consistent 
with Fig. 6 in for real- and redshift-space dark mat- 
ter and halo samples of similar bias if one compares the 
expected zero-crossing Tc position for their adopted cos- 
mology. The level of systematics should be small enough 
to allow accurate measurements of r^- We estimate that 
the error in ^(r) should be smaller than ^ 10"'^ (3 x 10"'') 
for 4(1)% error in Tc- 

Several groups have reported the absence of negative 
signal, and hence zero-crossing up to ~ 250 /i~^Mpc (see 
[2(J . l3ll | and references therein). It was speculated in 
dZ| that the large-scale SDSS £,{s) is affected by intrinsic 
errors and volume-dependent systematic effects. The im- 
plications of a possible constant systematic shift in the 
redshift-space clustering due to unknown observational 
systematic errors has been pointed out in [7,]. It was ar- 
gued in [l^l that lack of zero crossing in SDSS was simply 
due to cosmic variance, not unknown systematics. Inter- 
estingly, the WiggleZ redshift-space ^(s) at z = 0.6 Q, al- 
beit with large uncertainties, displays a crossover around 
~ 120/i~-'Mpc and negative signal above this scale up 
to 180 /i^^Mpc, as expected for current matter density 
observational constraints. 

We believe that the detection of the zero-crossing in 
the galaxy correlation function may offer some advan- 
tages as compared with finding the maximum in -P(fc), 
though mathematically both are equivalent. The indi- 
vidual Fourier modes in P{k) split the clustering signal 
in fragments and as a result the signal-to-noise and ac- 
curacy are expected to be smaller as compared with the 
cross-correlation signal around the zero-crossing which is 
coming from the convolution of the entire P{k) without 
any binning or splitting of the signal. 

To summarize, adding zero-crossing estimates will im- 
prove the accuracy of Wm because is more sensitive 
to Win than the BAO peak and because it does not re- 
quire the elaborate reconstruction methods used for BAO 
analysis [33^ ,34] . 
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